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1. Introduction. 



Standandard Hamiltonian mechanics in its homogeneous formulation is applied to the study of 
discontinuities representing rapid changes of Hamiltonians. Different formulations of Hamiltonian 
p^ ' mechanics are reviewed. An original representation of a positive homogeneous Hamiltonian system by 

^s) . an outer oriented coisotropic submanifold of the phase space is proposed. 

5^ ' The present note is a simplified and updated version of the article 

.^ . Hamiltonian dynamics v^rith discontinuities 

by 
Fiorella Barone, Ernesto A. Lacomba and Wlodzimierz M. Tulczyjew 



^ , 2. Formulations of Hamiltonian mechanics of autonomous systems. 

^^ , We present several equivalent formulations of Hamiltonian mechanics. In the analysis of discontinu- 

ities the different formulations play different roles. 
2.1. The traditional formulation. 



o 
^. 

^^ I In the simple traditional formulation of Hamiltonian mechanics the dynamics is generated by a 

smooth Hamiltonian 

H:P^R (1) 

^ |. defined on the phase space P. There is a symplectic form uj on P. The Hamiltonian generates a 

"^ '_ Hamiltonian vector field on P characterized by 

-(— > 

C^ ■ ixto^- dH. (2) 

Phase space trajectories of the system are integral curves of this field. They are solutions of the 
equation 

^; D^{p,eTP;p = X{Tp{p))} 

= [p, e JP ; V.peTP if ^Pi^P) = ^PiP)^ tl^cn {ujQ.pA Sp) = -{dH, Sp)} (3) 

generated by the Hamiltonian. 

If the configuration space of the mechanical system is an affine space Q modelled on a vector space V, 
then P is the affine space Q x V* . The Hamiltonian is derived by applying the Legendre transformation 
to a hyperregular Lagrangian 

L:Q xV ^R (4) 

and trajectories are solutions of the canonical equations 

. dHiq,p) . dHiq,p) 

2.2. The homogeneous formulation in terms of generating families. 

The homogeneous formulation is usually used with time dependent Hamiltonians. It seems to be 
indispensable when discontinuities ere considered. In the homogeneous formulation the phase space is 



replaced by the extended phase space P x K x K. Time and energy are included. The extended phase 
space is again a symplectic manifold with the symplectic form 

p*uj - de A dt, (6) 

where 

p:P xRxR^ P:{p,t,e)^ p, (7) 

t:P xRxR^R:{p,t,e)^t, (8) 

and 

e:PxRxR^R:{p,t,e)^e (9) 

are the canonical projections. We will use the afhne extended phase space P x R x M with P = Q x V* , 
projections 

q:P xRxR^V:{q,p,t,e) ^ q, (10) 

p:PxRxR^V*:(q,p,t,e)^p, (11) 

t:PxRxR^R:{q,p,t,e)^t, (12) 

and 

e:P xRxR^R:{q,p,t,e)^e, (13) 

and the symplectic form 

dp A dq - de A di. (14) 

The homogeneous Lagrangian 

L:Q xRxV xR+ -^ R: {q,t,q ,t') ^L L, ^ j t' (15) 

derived from a Lagrangian 

L:QxV^R (16) 

is not hyperregular. The Slow and Careful Legendre Transformation [1] must be used. The result of 
this Legendre transformation followed by a reduction operation is the family 

G:Px Rx]Rx]R+ -^ R: {q,p,t,e,T) ^ {H{q,p) - e) r (17) 

of functions defined on fibres of the fibration 

C.PxRxRxR+^ P xRxR: {q,p,t,e,T) ^ {q,p,t,e). (18) 

Extended phase space trajectories are solutions of the differential equation 

E^ l{q,p,t,e,q',p',t',e') ePxIRxMxV^xV^*xMxIR; t' > 0, e^ H{q,p), 

generated by the family (G, C) in the sense that {q,p, t, e, q' ,p' , t' , e') is in E if for some r > and all 
variations {5q, 6p, 6t, 6e, 5t) ^V xV* xMxlRxIR the equality 

{p',5q) - {Sp,q'} - e'St + t'5e = - (^^^o) - Up^ ^) +^'^^ " (^('^'^) " ""^ ^^ (20) 



holds. The equation E is positive homogeneous in the sense that if {q,p,t,e,q' ,p' ,t' ,e') G E, then 
{q,p, t, e, kq', kp' , kt' , fee') G E for each fc > 0. It foUows that parameterization of solutions is arbitrary 
except for orienation. If time is chosen as the parameter, then the equation assumes the form 

{q,p,t,e,q',p',t',e') ePxRxRxVxV*xRxR: t' = 1, e = H{q,p), 

equivalent to the equations (4) with the equalities 

e^H{Tp{p')), e' = (22) 

added. 

The family (G, C) is not the only family generating the dynamics. The factor H{q,p) — e in the 
definition (17) can be replaced by F{q,p^t,e), where F is any function on the extended phase space 
with two properties: 

(1) F vanishes on the set 

N ^ {{q,p,t,e) eP xRxM.; e = H{q,p)] (23) 

(2) the derivative 

dF 

^ '"" 

is positive at points of N . 
The result is a new function 

G:Px M xIRx]R+ ^ M: (q,p,i, e, r) i-^ F{q,p,t,e)T. (25) 

The family (G, C.) is again a generating object for the dynamics. 

2.3. Outer oriented coisotropic submanifolds. 

The set A^ is a coisotropic submanifold of the extended phase space. This set is practically the graph 
of the Hamiltonian. The characteristic distribution of N is the differential equation. 

E ^ l{q,p,t,e,q',p\t',e') ePxRxRxVxV*xRxR; e^ H{q,p), 

almost the same as the equation E. Only the orientation specified by i' > is absent. A method for 
restoring the orientation consists in specifying an outer orientation of N. Since iV is a hypersurface, 
an outer orientation is specified by assigning to each element {q,p,t,e) ^ N an equivalence class of 
vectors {6q, Sp, 5t, Se) G F x V* x M x IR transverse to TV. Within the equivalence class one can multiply 
vectors by positive numbers and add vectors tangent to A^. The correct orientation for our purpose 
is obtained by assigning to each element (g,p, i, e) of N the set of vectors Sq,Sp,St,5e satisfying the 

inequality 

dF dF dF OF 

(_.« + (,p_) + _« + _fc,>0, (27) 

where F is one of the functions introduced above. This set is the class of the vector {dq, 6p, St, 6e) = 
(0, 0, 0, 1). Let O be the outer orientation constructed in this way. The equation E is the set 

E ^ {{q,p,t,e,q',p',t',e') ePxRxRxVxV*x'RxR; {q,p,t,e) G N, and 

(p', 6q) - {5p, q') - e'St + t'Se > 

for each (Sq, Sp, St, Se) in the orientation O at (q,p, t,e)} . (28) 

The oriented coisotropic submanifold {N, O) provides an alternative to the generating families (G, C) • 
Uniquenes and simplicity arc the advantages of this representation of dynamics. 



3. Smooth models of discontinuities. 

We start with the smooth function 



,0 for s < 

exp(— s J tor s > 0. 



This function and its derivative 

for s < 

s^^exp(— s^^) for s > 

assign positive values to positive values of the argument. We observe that 



D(^:M^M:s^<{ _^ ^ ^ ^ (30) 



and 

It follows that 
The function 

has the values 



for s «; -1 

^^' + '^ = < exp(-(l + .)-)>0 for.>-l ^''^ 

exp(-(l-s)-i)>0 fors<l 

^(1 - ■«) = i n f -> 1 ^^^> 

tor s ;^ 1. 



ip{l + s) + ip{l - s) > 0. (33) 



ip{l + s) + ip{l - s) 



(34) 



X{s) = -1 for s s; -1, (35) 

X(s) = 1 for s ^ 1, (36) 

-1 < x(s) < 1 for -1< s < 1. (37) 



The derivative 



^ Dy(l + s) + D^(l - ,s) _ (^(1 + s)- y(l - s))iBif{l + s)- Dy(l - s)) 
^^ ' (p{l + s) + ip{l - s) {(p{l + s) + ^(1 - s))2 

_ (1 - x(g))Dy(l + -s) + (1 + x(g))D'^(l - s) 
tp{l + s) + ip{l - s) 



(38) 



is positive in the interval — 1 < s < 1. It follows that the funtion x is monotonically increasing in the 
interval — 1 < s < 1. The function 

/2s-2a\ /2a' -2s 

TCP w f 2s-a'-a \ ^ (^ ^3^ j " ^ 1^ a' - a , 

l"'" I \ a — a J Is — la \ / 2a — Is ^ 



a — a J \ a — a 

with a' > a, is a general purpose model of a step function. It has the values 

^ ^ f — 1 for s ^ a ^ ^ 

Va')(^) = |l forOa' ^'"^ 

and is monotonically increasing in the interval a < s < a' . We have 



^"^ X. Ja + k{a' - a)) ^ oo (41) 



for < A; < 1. 



4. A model of a simple Hamiltonian dicontinuity. 

A discontinuous Hamiltonian is an idealized representation of a smooth Hamiltonian rapidly changing 
across a hypersurface in the phase space. We will discuss a model of a simple discontinuity representing 
a monotonic rapid change of a Hamiltonian. 

Let Q be the Euclidean afhne configuration space of Newtonian mechanics of a single particle of 
mass m. The model space of Q is a vector space V with a Euclidean metric represented by a metric 
tensor g:V ^r V* . The phase space is the affine space P = Q x V* . We consider a hyperplane 

R={{q,p)eP; A{q,p)=0}, (42) 

where A is the function 

A:P^R: {q,p) ^{b,q- go) + {p - Po, «) (43) 

constructed with a phase ((jo,Po)j a- vector a G V^, and a covector b G V* such that 

{g{a),a) + {b,g-l{b)) = l. (44) 

Let H^: P ^M. and iJ+: P ^ M be smooth Hamiltonians. The function 

Hs:P^^:{q,p)^\{l-X^_,^S)^A{q,p)))H-{q + g-\b)5,p + g{a)5) 

+ ^ (l + X(-^,.-)(^(«'P))) ^^('^ - 9-\b)5,P- 9{a)5) (45) 

is a model of a discontinuous Hamiltonian designed for studying trajectories of the particle originating 
in P~ and directed towards R. 

In the extended phase space the dynamics associated with the model Hamiltonian is the differential 
equation 

Es = |((7,p,i,e,g',p',t',e') gPxMxKxVxV^xKxM; f'>0, e = Hs{q,p), 

generated by the family {Gs , C) composed of the function 

G^tPxM xM xM+ ^K:(g,p,i,e,T) ^ {Hs{q,p) - e) r (47) 

and the fibration (. introduced earlier. The set Eg is the oriented, or directed, characteristic distribution 
of the outer oriented coisotropic submanifold {Ns, Os) with 

Ns = {{q,p,t,e)ePxRxR: e^Hs{q,p)} (48) 

and the orientation Os determined by F{q,p,t,e) ~ Hs{q,p) — e. 

5. The limit of the model. 

In the limit oi 6 = the model Hamiltonian Hg defines a singular Hamiltonian 

H:P\R->R (49) 

equal to H~ on 

P- ^{{q.p)(.P; A{q,p)<Q}, (50) 
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equal to H^ on 

P+ = {iq,p)eP; Aiq,p)>0}, (51) 

and not defined on the hyperplane R. The traditional formulation of Hamiltonian mechanics can 
be applied to the model Hamiltonian Hg and limits of trajectories can be studied. The traditional 
formulation can not be applied directly to the singular Hamiltonian H. This is why the homogeneous 
formulation should be used. 

The limit of Ng is the union of coisotropic submanifolds 

A^" = {{q,p,t,e) gPxMxK; {q,p) E P' ,e = H' {q,p)} , (52) 

N+ = {iq,p,t,e)ePxRxR; {q,p) E P+ ,e = H+ {q,p)] , (53) 

and 

M ^ {{q,p,t,e) e P xRxR; {q,p) e R, 

H^{q,p) < e < H+{q,p) or H+{q,p) < e < H^{q,p)] , (54) 

and the two edges 

M" = {(g,p,t,e) e P xRxR; {q,p) e R,e = H' (q,p)} (55) 

and 

M+ ^ {{q,p,t,e)eP x^x^- {q,p) e R,e = H+{q,p)} . (56) 

The union 

N ^ N- U N+ U M U IVr U M+ (57) 

is the limit of Ng in the sense that for each neighbourhood U of {q,p, t,e) (z N the intersection NsOU 
is not empty for some 5 > 0. Outer orientations of the coisotropic components N^, N^, and M 
are derived from the outer orientation Os of Ng. Outer orientations O^ and 0+ of iV~ and N'^ are 
determined by the functions H~{q,p) — e and H'^{q,p) — e respectively. The orientation O of M at 
{q,p,t,e) is determined by the function A if H^{q,p) > H^{q,p) and by —A if H^{q,p) < H~{q,p). 
Characteristics of {N^ , O^) and (A^+, 0+) are easily constructed from solutions of Hamilton equations 
for H~ and H^. Let an oriented characteristic of {M,0) start at the point (qo,po,to,H^{qo,po)) E 
M~. The curve 

s ^ (qQ + sa,po - sb,to,H^{qo,po)) (58) 

is a suitable parameterization of the characteristic if H^(qQ,pQ) > H^ {qQ,pQ) and the curve 

s ^ {qo - sa,po + sb,to,H^{qo,po)) (59) 

is a suitable parameterization of the characteristic if H^(qo,po) < H^ (qo,po)- The parameter is 
restricted by the inequalities in the definition (54) of M. 

We analyse the limit of a trajectory of the particle passing through the region 

Rs^{iq,p)eP: -S<Aiq,p)<S}, (60) 

where the discontinuity develops. We consider a trajectory originating in 

Pg- ^ {iq,p) e P ; A{q,p) < -6} (61) 

and entering the region Rs at (go — fJ^^{b)S,po — g{a)S) at t ~ to with energy 

e = Hsiqo - g-\b)d,po - gia)d) = H-{qo,po). (62) 
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We assume that the condition 

, dH- ^ A / dH- ^ , \ 

b,-g:^{qo,Po)) -{-g^(qo,Po),a) >0 (63) 

is satisfied. This condition excludes from the analysis trajectories tangent to the hyperplane R at 
the point of inpact. The particle will encounter no discontinuity unless H^{q,p) ^ H^{q,p) in a 
neighbourhood of the point of inpact. In the homogeneous formulation trajectory appears as oriented 
characteristic of the oriented coisotropic subnianifold {Ns,Os). The continuation of the characteristic 
in the region Rg is well represented by the time parameterized curve 

t^{js{t),t,H-{qo,po)), (64) 

where 

js-.'R^Rs (65) 

is an integral curve of the Hamiltonian vector field 



V r> 1^ T^* / \ f^^S dHs 



(66) 



Derivatives of the model Hamiltonian Hs are the functions 

dHs ^ ,.* . ^ 1 /, ... x^\ dH- 



dq 



:Rs^V*:{q,p)^-[l-X^,^,^{A{q,p)))-^{q + g-'{b)S,p + g{a)5) 



+ 2V + ^(o,.)(^(«'P))) -Q-^l- a-\mP- 9{a)S) 



dq 



+ i (H+iq - g-\b)5,p- g{a)5) -H-{q + g-\b)5,p + g{a)S)) Bx^^,^iA{q,p))b (67) 



and 



^■.Rs^V:{q,p)^l{l-X^,^sM{q,p)))^{q + 9-Hb}S,p + g{a)S) 

m_ 

dp 
+ i {H+{q - g-\b)5,p- g{a)S) - R- {q + g-'(b)5,p + g{a)S)) Dx^^ ,)(A(g,p))a (68) 



1 / \ rlH+ 

+ - (l + X(o,,)(A(<z,p))) -g-{q- 9-'{b)S,p- g{a)5) 



The factor 

K{q,p) = i {H+{q - g-\b)6,p - g{a)5) -H-{q + g-\b)S,p + g{a)S)) Bx^^^^^{A{q,p)) (69) 

diverges to plus or minus infinity at S — 0. The curve (64) is an integral curve of the vector field 

Xs- Rs -x ^ >^ R ^ V X V* X RxR:{q,p,t,e) ^ (^5,1,0). (70) 

An equivalent representation of the characteristic is obtained by replacing the curve (64) by an integral 
curve of the renormalized field 

Zs:RsxRxR^V xV* xR xR: {q,p,t,e) ^Xs{q,p,t,e)\\K{q,p)\\-'^. (71) 

The field 

Z:RxRxR^V xV* xR xR: {q,p,t,e) ^ K{q,p)\\K{q,p)\\-\a,-b,0,0) (72) 
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is the limit of Z5 . The curve 

s 1-^ (go + sK{qQ,pQ)\\K{qo,pa)\\^'^a,po - K{qo,po)\\K{qQ,pQ)\\^'^b,tQ, H^ {qo,po)) (73) 

is an integral curve of Z starting at {qo,po,to, H^ {qo,po). It is a representation of the limit of the 
continuation of the analysed characteristic as a characteristic of {M, O). 

Ineqalities in the definition of M restrict the values of the parameter s to an interval ]0, si[. The 
final point on the curve (73) corresponding to s = si is either a point {qi,pi,ti, H^ {qi,pi) in M^ or 
a point (gijpi, ti, H^{qi,pi) in M"*". Equalities ti = to and H^{qi,pi) = H^{qi,pi) = H^{qo,po) are 
always satisfied. If the final point is in M~ and the condition 

^'^(9i'Pi))-(^(9i,Pi),«y <0 (74) 

is satisfied, then the particle is reflected back into P~ and continues as a characteristic of {N^ , O^) 
with {qi,pi,ti, H~{qi,pi)) as the initial point. If the final point is in M"*" and the condition 

^^^('Zi,Pi)y-^^^(9i,Pi),«y >0 (75) 

is satisfied, then the particle enters P+ and continues as a characteristic of {N^, 0+) with 
(qi ,pi,ti, iJ+ (qi ,pi)) as the initial point. The continuous passage between the characteristic of (M, O) 
and a characteristic of {N" , O^) or (N^, 0+) with only one point added is proved as it was done above 
for the initial point {qo,po,to, H^{qo,po)) of the characteristic of {M,0). 

6. A Transition Principle. 

The following principle summarizes the results of the preceding section. It is formulated in terms of 
definitions introduced in this section. Transversality conditions (63) and (74) or (75) are assumed to 
hold. 

The Transition Principle. The Hamiltonian discontinuity reflecting a rapid monotonic change of a 
Hamiltonian is represented by the union 

N ^ N~ U N+ U Ad U M- U M+ (76) 

of coisotropic submanifolds N^ , N^ , and M, and the borders M^ and M+ of M . The coisotropic 
components of N have their natural exterior orientations. The trajectory of a particle crossig the 
discontinuity is a continuous arc in the extended phase space composed of directed characteristics of 
the oriented coisotropic components of N completed with isolated points of the borders. 

7. A simple example. 

We consider a Hamiltonian with a discontinuity due to a monotonic rapid change of the potential 
from a constant value U^ to a constant value U^. There is no discontinuity unless U^ ^ U^ . The 
discontinuity occures on the hypcrplane 

R = {qeQ; Biq)=0}xV*, (77) 

where B is the function 

B:P ^R:q^ {b,q-qo) (78) 

constructed with a configuration qo and a covector b G V* such that 

{b,g-lib))^l. (79) 



The function 



Us:Q^R:q^^{l- X(„ .^ {B{q))) [/- + 1 (l + x^„,^ {B{q))) U+ (80) 



is a model of the discontinuous potentiaL The function 

Hs: P^R: {q,p) ^ ^{p.g-^p)) + Us{q) (81) 

is the corresponding model of a discontinuous potential. 
The system is represented by the union 

iV = iV"UA^+UM (82) 

of coisotropic submanifolds 

N- = Uq,p,t,e) ePxRxR; {q,p) eP-,e = l-{p^g-\p)) + U-\ , (83) 

N+ - |(g,p,t,e) e P X M X M ; {q,p) EP+,e= l-{p, g-\p)) + U+\ , (84) 

and the closure 



M ^ \ {q,p,t,e) e P xRxR: {q,p) e R, 
Zm Im 



of the coisotropic submanifold 



or^{p,g-Hp)) + U+ ^ei^^{p,g-\p)) + U-\ (85) 

Zra Zm J 



M= <^ {q,p,t,e) e P xRxR; {q,p) e R, 

7^{p,9-Hp)) + U- <e< -^(p,5-i(p)) +[/+ 
Zm Zm 

OT ^{p,g-\p)) + U+ < e < ^{p,g-\p)) + U-\ (86) 

The coisotropic components arc oriented according to the rules stated in Section 2.3. A directed 
characteristic starts in P^ x M x R and proceeds in the direction of the discontinuity. The time 
parameterized description is given by 

t^ lqo + -g'Hpo)t,Po,t,^{po,g-^{po)) + U-) (87) 

\ m Zm / 

with t < 0. The condition {b, g^^{po)) > is satisfied. The characteristic meets the discontinuity at 

qa,Pa,0,—{po,g^^{po)) + U' 
The characteristic continues as a directed characteristic of (M, O). If C/+ > U^ , then 

qo,po- sb,0,^{po,g-\po)) + U-j (89) 
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is a convenient parameterization of the characteristic. The value of the parameter s is restricted to 
the interval in which the inequality 

0^-{b,g-\po))s-^s^^U+-U- (90) 

m Zm 



holds. The graph of the polynomial 



S{s)^-{b,g-'{po))s~^s' (91) 

m Zm 

is a parabola and the maximum value of S{s) is 

^ib,g-\p„)r (92) 

at 

■s = {b,g-'{po)). (93) 

The parameter of the initial point is the root s — so = oi the polynomial. We will find the parameter 
s = si of the terminal point. We consider three cases. 



(1) If 



T^{b,g-\po)f<U+-U-, (94) 

Zm 

then the boundary M+ can not be reached. The terminal point 

(qo,Pi,0,^ipi,g-\pi))+uA (95) 

with 

Pi=Po-2{b,g-\po))b (96) 

is in M" and its parameter is 

s^si^2{b,g-'{po)). (97) 

This is the other root of the quadratic polynomial S{s). 

The component of the momentum orthogonal to b is not changed in the collision since 

Pi - (6, g-Hpi))b = Po- {b, g-\po))b. (98) 

The component parallel to b inverts its direction: 

{b,g-\p2))b = -{b,g-\p,))b. (99) 

After the collision the trajectory is again in P^ x R x M and is described by 

t^ (qo + -g'\pi)t,Pi,t,:^{pi,g-\pi)) + U-) (100) 

\ m Zm J 

with t > 0. The point moves in P^ x M x M away from the discontinuity since 

-{b,g-\pi)) = --{b,g-\po)) < 0. (101) 

m, m 
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(2) If 



^{b,g-\po)f>U+~U-, (102) 



then the terminal point is in M+ and the trajectory continues across the discontinuity into 
P+ X R X M. The parameter s = si of the terminal point is the smaller of the two solutions of 
the quadratic equation S{s) = C/+ — U^ . The value of this root is 



si = {b,g-\po)) - V{b,9-HPo))^ - 2m{U+ - C/-). (103) 

The extended phase 

(qo,Pi,0,^{pi,g-Hpi)) + uA (104) 

with 

Pi=Po- s,b = po- (^{b,g-\po)) - V{b,g-'{Po))^ - 2miU+ - U-) ) b (105) 

is the terminal point. The component of the momentum orthogonal to b is again conserved. 
The component in the direction of b is reduced from 

{b,g-\p^))b (106) 

to 

{b,g-\pi)) b = ^/{b,g-^PoW - 2m{U+ - U-) b. (107) 

The energy is conserved. The characteristic continues into P+ x M x M and is described by 

t^ {qa + -9-\pi)t.Put,^{p^,g-\p^)) + U+] (108) 

with i > 0. The point moves in P+ x M x M away from the discontinuity since 

-{b,g-\pi)) = -^{b,g-^po))^-2miU+-U-) > 0. (109) 

m m 

(3) The case 

^{b,g-\po)? = U+-U- (110) 

must be considered as approximated by one of the two cases above. If it is approximated by 
case (1), then the terminal point of the characteristic arc is in M^ and its parameter is 

s = si = 2(6,5-i(po)). (Ill) 

The outcome is a reflection of the trajectory back into P^ x M x M as in case (1) above. If this 
case is approximated by case (2), then the terminal point is in M^ and its parameter is 

5 = si-(6,g-i(po)). (112) 

We do not analyse this possibility since the transvcrsality requirement (75) is violated. 
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We assumed initially that t/+ > U^ . Wc now proceed with the analysis of the case C/+ < U^ . The 
initial trajectory starts in P^ x K x M and has the form (87). The characteristic arc of (M, O) starting 
at (q,p,t,e) = {qo,po,0,eo) with s = sq = is described by 

s^ (qo,Po + sb,0,^{po,g-\po)) + uA (113) 

The value of the parameter s is restricted to the interval in which the inequality 

^ L{b, g-\po))s + ^s^ ^ U- - U+ (114) 

m Zm 

holds. The terminal point of the characteristic arc corresponds to the value of the parameter obtained 
as the positive solution 



si = V{b,9-'iPo))^ + 2m{U- - U+) ~ {b,g-\p^)). (115) 

of the quadratic equation 



The extended phase 



with 



-{b,g-\po))s + ^s^ = U- -U+. (116) 

TO Im 



qo,Pi,0,^{pi,g-^(pi)) + U+] (117) 

2m ' 



Pi=Po + sib = po + (^^{b,g-^{poW + 2m{U--U+) - {b, g-\po)) ) b (118) 

is the terminal point. 

The component of the momentum orthogonal to b is conserved. The component in the direction of 
b is increased. The energy is conserved. The trajectory continues into P+ x M x R and is described by 

t^ (qo + -g-\pi)t,pi,t,^{pi,g-\pi)) + U+) (119) 

with t > 0. The point moves in i?+ x M away from the discontinuity since 
1 ,. _i, ., 1 



{b,9-\pi)) = -V{b,9-HPo)y+2m{U--U+) > 0. (120) 

in m 

8. Non monotonically changing Hamiltonians. 

A generalization of the Transition Principle is needed to handle discontinuities of Hamiltonians 
reflecting non monotonic rapid changes. We give a simple examle of such discontinuity. Let t/", U^, 
and C/^ be constants such that U^ > Uq and U^ > U2- Let B be the function introduced in the 
preceding section and let 

<5o = 0, h^^S, 52^5. (121) 

The function 

Us-Q ^M.:q^ < ^ , (122) 

is a non monotonically changing model of a discontinuous potential and 

Hs: P ^ R: {q,p) ^ ^{p,g-\p)) + Us{q) (123) 

Zm 
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is the corresponding model Hamiltonian. If a particle starts in P x M x M and follows the time 
parameterized characteristic 

t^ (qo + -9'HPo)t,Po,t,^{po,g-\po)) + uA (124) 

\ m Zm J 

with t ^ and {b, g^^{po)) > then two outcomes of the collision at 

(qo,Po,0,:^{po,9-Hpo)) + uA (125) 

can be considered: 

(1) If 

-L{b,g-\po)f<U'-U", (126) 

Zm 

then the continuation of the characteristic is represented by 

s^ (qo,Po-sb,0,^{po,9-\po)) + uA (127) 

with ^ s ^ 2{b, g^^{po)). This characteristic arc terminates at 

(qa,Pi,0,^{po,9-\po)) + uA (128) 

with 

Pi=Po-2{b,g-\po))b. (129) 



The particle reenters the region P x R x M and follows the characteristic 

t^ (qo + -9-Hpi)t,Pu 7^{put,g^\pi)) + uA 
\ m Zm J 



(130) 



with i > 0. 
(2) If 

^{b,g-\po)?>U'-U\ (131) 

Zm 

then the trajectory of the particle in the extended phase space has four components. There is 
the original trajectory (124) followed by the characteristic arc 

s^ (qo,Po-sb,0,^{po,g-\po)) + uA (132) 



with =^ s ^ {b,g^^{po)) - \/{b,g^^{poW - '2m{U'^ - U°)- This arc terminates at 

'qo,Pi,0,^(po,g-\po)) + uA (133) 



with 

pi^po- {b,g-\po))b+^{b,g-Hpo))^-2m{m - U^) b (134) 

The next component is the arc 

s^ (qo,Pi+sb,0,^{pi,g-\pi)) + U'^ (135) 

13 



with < s ^ A/(6^5^^T(pI)pT"2m(C7^~~C72y - {b,g-^{pi)). This arc terminates at 

[ciu,P2,0,^(j>2,g-\v2)) + U^\ (136) 

with 

P2=Pi'- {b, g-\pi))b + ^{b,g-^{p,))^-2m{U^-m) b 



= Pa-{b,g-'ipo))b+^{b,g-^ipo)y-2miU^-m)b. (137) 

The particle enters the region P+ x M x M. The last component is the characteristic 

t^ (qo + -g'HP2)t,P2,^{p2,t,g-Hp2)) + uA (138) 

\ m 2m J 

with t > 0. 

The example shows that the Transition Principle formulated in Section 6 is incomplete. A suitably 
extendet principle must be formulated if more complex rapid changes of the Hamiltonian are to be 
taken into account. 

9. The Transition Principle. 

Vinogradov [2] [3] claims that Hamiltonian mechanics does not provide information on how a tr- 
jectory of a particle is affected by the presence of a dicontinuity in the Hamiltonian. The Transition 
Principle postulated by Vinogradov is expected to fill this gap. We are using simple examples to reveal 
the defects of this Principle. 

We consider a symplectic phase space P composed of two open components P^ and P+ separated 
by a hypersurface R. There are two smooth Hamiltonians H^ and H^ defined on P. We introduce 
the closed sets P — P~ U R and P = P+ U R corresponding to $_ and $+ in [4]. A discontinuous 
Hamiltonian is a function H on P^ U P+ constructed by setting H\P~ = H~ and H\P^ = iJ+. The 
Hamiltonian H is not defined on R. The Hamiltonian vector fields generated by the Hamiltonians H~ 
and 7?"*" will be denoted by X~ and X^ respectively. The hypersurface _R is a coisotropic submanifold 
of the symplecic phase space. It is foliated by one-dimensional characteristics. The characteristic 
passing through a point p E R will be denoted by jp. 

We will consider the trajectory of a particle originating in P~. We extract from [4] definitions of 
concepts applicable to this situation. 

Definition 1. A point p G P is called an in-point for iJ+ if X^(p) is directed towards P . A point 
p G P is called an in-point for H^ if X^{p) is directed towards P . 

Definition 2. A point p' G jp is said to be decisive for {p, H^) if it is either an in-point for H^ and 
H^{p') = H^{p) or it is an in-point for H^ and H^{p') = H^{p). 

The following formulation is the one found in [4] cxept for the symbols used. 
The Transition Principle. When the moving phase point reaches the hypersurface R from P at a 
point p, its trajectory is to be prolonged, starting from any decisive point for {p,H^), as the trajectory 
of the corresponding Hamiltonian. 

The following formulation is more explicit. 
The Transition Principle. When the moving phase point reaches the hypersurface R from P at a 
point p, its trajectory is to be prolonged, starting from any decisive point for {p,H^). If the decisive 
point is an in-point for H^ , then the trajectory is to be prolonged in P as the trajectory of H^ . If 
the decisive point is an in-point for II~ , then the trajectory is to be prolonged in P as the trajectory 
ofH-. 
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I will interpret the phrase "is to be prolonged, starting from any decisive point" as meaning that a 
prolongation trajectory is to start at each decisive point and not just one of the decisive points. This 
is explicitly stated in [3]. 

We will apply the Transition Principle to the simple example of discontinuity of Section 7. The 
nature of the discontinuity as reflecting a rapid change (monotonic or not) of the potential is not 
specified. The analysis applies to the example in Section 7 with U" = C/" and C/+ = U"^. The two 
Hamiltonians H~ and H'^ generate the same vector field 



X:QxV*^QxV*xVxV*: {q,p) ^ {q,p,q,p) = (q,P, -g"'(p),o') 

\ m J 



(139) 



The parameterized characteristic through (<7,|?) G i? is the curve 

7(,,p):IR^i?:s^(g,p + s6). (140) 

A particle starts in P^ and follows the trajectory 

t^ (go + -5"'(Po)t,Po) (141) 

\ TO / 

with {b,g~^{po)) > and t ^ 0. At i = the particle enters in contact with the discontinuity at 
{QojPo)- The trajectory will be prolonged for t ^ starting from each decisive point for {{qo,Po), H~). 
A point {qi,pi) — {qo,Po + sb) is decisive if 

:^{pu9-\pi))^:^{po,9-HPo)) (142) 

and 

-(6,5"i(pi))<0 (143) 

m 

or 

-^{pi,g-\pi)) + ^U = ^(j>^,g-\p^)) (144) 

Zra Zm 

and 

-{b,g-\pi))^0. (145) 

TO, 

In the first case the new trajectory reenters P . In the second case the trajectory continues into P . 
A point {qi,pi) = {qa,Po + sb) will satisfy the equation (142) if the parameter s is a solution of the 
quadratic equation 

s^ + 2s{b,g-\po))^0. (146) 

There are two solutions 

s^O and s = -2{b,g-\pa)). (147) 

Of the two corresponding points 

(<Zo,Po) and {qo,po-2{b,g-\pi))b) (148) 

only the second satisfies the condition (143) and is decisive. The new trajectory 

q = qo + -g-\po-'2{b,g-\po))b)t, p ^ po - 2{b,g-\p„))b (149) 

TO, 

with t ^ reenters P . Note that the reflected trajectory is predicted even if AU = 0. In the absence 
of a discontinuity the choice of the hyperplane R is arbitrary. Thus the reflections predicted by the 
Transition Principle occur at each point along the trajectory of a particle and in all directions. 
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A point {qi,pi) — {qo,Po + sb) will satisfy the equation (144) if the parameter s is a solution of the 
equation 

s'^ + 2s{b,g-\po))+2mAU = 0. (150) 



There are no solutions if 



l-{b,g-^{po)f<AU (151) 



since the discriminant 

^{{b,g-'{po)f-2mAU) (152) 

is negative. There is one solution 

s = -{b,g-\po)) (153) 

if 

^{b,g-\p,)f = AU. (154) 

The corresponding point 

{qo,Po-{b,g-Hpo))b) (155) 

is decisive. The vector field X is tangent to R at this point. The inequality (145) is satisfied. The new 
trajectory 

g = go H g^^ipo- {b,9^^{Po))b)t, P = Pa- {b,9^^{po))b (156) 

m 

with t J^ is in i?. There arc two solutions 

s = -{b, g-\pa)) - V{b,9-HPo))^-2mAU (157) 

and 



s = -{b,g-\p„)) + V{b,9-'iPo))^ - 2mAt/ (158) 

if 

^{b,g-\po)f>AU. (159) 

Of the two corresponding points 

(go,Po - {b,g-\po))b - V{b,9-^{PoW - ^mAU &) (160) 

and 

(qo,Po - {b,g-\po))b+V{b,9-'{Po))^ - ^mAU fe) (161) 

only the second satisfies the condition (145) and is decisive. The new trajectory 

q^qo + -g-^ (po - {b,g-\pa))b + ^{b,g-\pa))^ -2mAU b) t, 

p = Po-{b,9-\pa))b+^{b,g-^po))^-2mAUb (162) 

with i ^ continues in P+. A split trajectory is obtained since there is already a reflected trajectory. 
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10. Concluding remarks. 

We have applied tlic homogeneous formulation of standard Hamiltonian theory to the analysis of 
discontinuities treated as representing rapid changes of Hamiltonians. We find that if some information 
about origin of a discontinuity is provided, then the behaviour of a particle encountering the discon- 
tinuity is fully predicted by the standard Hamiltonian theory without external postulates. We find 
that, in its present formulation, the Transition Principle postulated by Vinogradov is not valid since 
it predicts split trajectories and un-provoked reflections. In a modified version the Principle could be 
used as a definition of elastic transitions across discontinuities. 

We propose the following version of the Principle: 

The Transition Principle. When the moving phase point reaches the hypersurface R from P at a 
point p, its trajectory is prolonged, starting from one of the decisive points for {p, H^). If the decisive 
point is an in-point for H^ , then the trajectory is prolonged in P as the trajectory of H^ . If the 
decisive point is an in-point for H^ , then the trajectory is prolonged in P as the trajectory of H^ . 
We propose the following definition: 

Definition 3. The transition of a trajectory across a discontinuity is non dissipative, or elastic, if 
the above Transition Principle holds. 

Some restrictions were used in our analysis: we treated only time independent Hamiltonians with 
time independent discontinuities, we assumed affine geometry for the phase space, we imposed transver- 
sality conditions (63), and (74) or (75). The first two restrictions are easily removed. An analysis of 
trajectories not satisfying the transversality conditions is still to be done. 
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